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We analyze the decays B° ctj^^ B^'" — > f(,K~'^ and show that within the factorization 
approximation a phenomenological consistent picture can be obtained. We show that in this ap- 
proach the Oe operator provides the dominant contributions to the suppressed channel — » ajvr". 
When the ao is considered a two quark state, evaluation of the annihilation form factor using Per- 
turbative QCD implies that this contribution is not negligible, and furthermore it can interfere 
constructively or destructively with other penguin contributions. As a consequence of this ambi- 
guity, the positive identification of — > i^'^a.Q can not distinguish between the two or four quark 
assignment of the ao 0]. According to our calculation, a best candidate to distinguish the nature of 
ao scalar is Br{B~ t^^o-q) since the predictions for a four quark model is one order of magnitude 
smaller than for the two quark assignment. When the scalars are seen as two quarks states, simple 
theoretical assumptions based on SU{2) isospin symmetry provide relations between different B 
decays involving one scalar and one pseudoscalar meson. 



INTRODUCTION 



B factories provides large samples of B — B mesons allowing the study of physicalphenomena such as CP violation, 
the determination of the CKM mixing angles and the search for new physics 0, Clearly, hadronic physics will 
beneficiate of the high statistics achieved, and the study of processes with small branching ratios will be possible. The 
full understanding of the B physics is still lacking as well as a systematic first principles description of the phenomena 
involved. Instead different theoretical approaches are compared to data and assumptions such as factorization, or 
estimation of the relative size of different contributions (tree level, annihilation, penguins, final state interactions) 
can be tested. This can be achieved in processes where the dominant contributions are suppressed by symmetry or 
accidental cancellations. 



The BABAR and Belle collaborations already reported precise measurements of non-leptonic B meson decays 
involving scalar mesons with branching ratio of order as low as 10^®. Thus for example, for the foK^ channel, 

besides the branching ratio the CP violating asymmetries are reported and, from the two pion spectrum, the authors 
are able to obtain the mass and width of the /o (980)0. This is not the case for the B ao(980)7r where branching 
fraction for given final states are reported -in particular Cq tt"*"- however in this case it is not possible to separate the 
from the & decays, unless a dominant decay mechanism is assumed^. In this context it is worth remarking 
that the B — > ao(980)7r decay was suggested as a place where a, the weak mixing angle, could be measured through 
the CP violating asymmetries 0|. However, it was shown that the B Oo^tt^ is suppressed by G parity and also 
by isospin, which implies that in the symmetry limit no CP violating asymmetry is expected to be experimentally 
accessible 0. Thus, theoretical arguments support the idea that the B^ o,qT^~ is strongly suppressed, so that 
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the reported branching ratio can be identified with the dominant B° ^-cT'""^ decay. As a by product, in |g the 
author concludes that the positive identification of / '^o^'"'^ is an evidence against the four quark assignment 
of ao or else, for the breakdown of perturbative QCD. Both, the suppression of the 5° — ^ o.q''^^ channel as well as 
the potential evidence in favor or against a four quark state are rather interesting observation that deserves further 
analysis. 

The low lying scalar sector of QCD represents a major challenge [^. From the experimental point of view, the 
determination of the nature of existing states has not been achieved while from the theory side no consistent interpre- 
tation of the experimental data exists 0. This is so even though a number of processes involving the appropriated 
final state in the kinematical region of interest have been analyzed. Thus, for example, <j) — > 77777, J ~* (/)7r7r, cjyKK 
and central production involve the /o(980) and ao(980) JJiJ, whereas the di-pion in the TjnS) — > T{mS)TnT and 
D TTTTTT decays include the kinematical region where the /o(600) is expected to appear [ul • Unfortunately, 

although different processes are included in the analysis, data is not good enough to provide a clear picture of the 
scalars. In fact no consensus exist even on the fundamental intrinsic properties (mass and width) of the low lying 
scalar mesons 0, • 

The appropriated theoretical description of non-leptonic B decays involving scalar mesons is important not only to 
understand the nature of the scalars but also because these must be considered as background to other processes of 
interest in B physics. Since scalars, vectors and tensors couple to two pseudoscalars, the following decays lead to the 
same final state: B — > PV, B SP, B — > TP and B PPP, where S, V, T and P stand for scalar, vector, spin 2 
and pseudoscalar meson respectively. B decays involving scalar mesons have been considered by a number of authors. 
Thus for example in |l4l | the tree level Hamiltonian and quark model calculation are used to predict branching ratios, 
while sum rules [isj and gluon-penguin dominance (6 —> sg) are the basis to interpret the scalars produced in B 
decays in terms of glue balls 0|, or still using QCD corrected Hamiltonian plus factorization to propose evidence 
for two quark content of the /o(980). 

In this note we analyze the fo{980)K and B ao(980)7r decays using the factorization approximation. To 

this end we use the AB = 1 weak Hamiltonian including QCD corrections to next to leading order. To evaluate the 
matrix elements we use values reported in the literature or model dependent estimates of the decay constants and 
form factors. In particular, the annihilation contribution is considered and evaluated using Perturbative QCD. This 
is important in estimating contributions previously neglected, and it is also relevant to quantify the statement in 
regarding the four quark nature of the ao(980). 



B° a±7rT AND 3°'' foK°'- 

FoUowing the conventional approach m El HO, m 113, we start with the AB — 1 effective Hamiltonian H^f f{q — 



d,s) 



K 



Gf 
V2 



h.c. 



(1) 



where Xq'q = Vq'bV*,y,q ~ d,s,q' — u,c,t,Vij are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The 
Wilson coefficients, Ci, including next to leading order QCD corrections, are evaluated at the renormalization scale 
H ~ ms/2. We use the conventions and Ci values reported in |2Cl| . It remains to evaluate the matrix elements between 
the states of interest. 



A{B'^PS)= <PS\Heff\B> 



(2) 



P and S stand for pseudoscalar and scalar meson respectively. In terms of the amplitude the branching ratio is given 
by: 



Br{B PS) ~ TB 



G%\A{B^PS)\' 
327r rriB 



(3) 
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with tb the appropriated B meson hfetime {Tg+ = 1.65 • 10 s, Tgo = 1.56 • 10 s). 

The matrix elements are evaluated using the assumption of factorization. In that approximation the matrix elements 
of interest are given by: 
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2{ae + as)X ° - 04 - -(og - ay) - -oio + 



(ag + a8)m|^s- 



mu{mb + md) 



+ a4 + aio + 



fh{mi, + m„) °-o 



(6) 



^B-^TT-afl - A„</ai(^B-aO +^^7r-) ~ Afd 



(ag + a8)m| 



+ a4 + aio - 



rh{mb + ffi) 



(.04 + aio - ^/ , ^ A-'^s-ag 
m[mi, + rriu ) " 



Xf„;_+(a8-2ag)X^l^_ 



(7) 



Ar-_ 



Ausftl 

+ a4 + aio - 



A(s 

2(ao + as)mj^ 



a-i + aio 



(m„ + ms){mb + m„) 



2(ag + a8)m| 
(rrifc + mu){ms + iriu) 



(2a6 + a8)X^=- 



fOK- 



(8) 



BO-*f°KO 



(a4 



flio 



(2ag - a8)m| 



a4 ■ 



aio 
2 



(to 



(mb + md){ms + md) 
(2ag - a8)m^ \ ^0 
,+TO.)(m, + TO.)j ^^'V" - (^"^ - 



(9) 



where to = (to„ + md)/2. For future reference in Table 1 we quote the numerical values of the a^ coefficients. These 
expressions are obtained by inserting the vacuum between the currents in all possible ways, and a typical product 
of matrix elements is parameterized in terms of form factors and decay constant as: 



X 



i:^^ = < Ky"\{sd)L\0 X 0\{db)L\B" > = fsimjo - mi,)F[ ^\m's) 



(10) 



In the appendix we define in detail all of the ^. Let us start by summarizing our knowledge about the decay 
constants and form factors entering the calculation. We can classify these in four categories: Pseudoscalar decay 
constants ( f^r, fw, /«). The values of the two former decay constants are taken from [4] while for the later we use 
/b = 170 MeV|23. The second kind are the scalar decay constants (/a,//)- For these we use pubhshed values 
estimated using theoretical arguments 0- We then have form factors of the type FQ^{m\), evaluated at the tob 
scale, i.e. calculable with perturbative methods. In Table l(Tl)) . we quote the values we use. Finally we need the 
form factor ^ [m^p) , ^ {m^) where S and P stand for a scalar or pseudoscalar meson. For the decay under 
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TABLE I: Numerical values for the effective coefBcientes a1 for 6 — > d transitions at scale ^ ~ mj, (for 03, aio in units of 
10-* (a2i-i = C2i-i + C2i/iV, a2i = C2i + Ca.-i/iV)) 
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TABLE II: Numerical values of the form factors 



consideration these are evaluated at most at the scalar scale (around 1 GeV), i.e. these can not be computed using 
perturbative methods and, in general, few is known about their values. 

For the B foK", B- f^R- and B cott we require F^'' , F^^ , Fq^\ F^'^l^F^'' and Ff'^. The two first 
(FfP'^ , Fq^) are relatively well-known and we shall use the value presented in refs.|2J|. In principle, the other four 
form factors could be determined using experimental data^ Q II] given in Table but since the experimental 
data involve large error bars, we prefer to evaluate some of these form factors. In this paper we focus our interest on 
the annihilation effects, for that reason we present an estimate of F^""^ using perturbative QCD and considering oq 
as a two or four quark state. 



ANNIHILATION FORM FACTORS FROM PERTURBATIVE QCD. 



At tree level the B^ strongly suppressed due to the absence of second class currents. To get an evaluation 

for such a decay an estimate of the contribution of the B annihilation is necessary. The annihilation amplitude is 

proportional to ^1^°-^+^ which itself is proportional to the Fq° ^ (?7t-^) form factor. At the scale m^, perturbative 
QCD provides an adequate framework to evaluate this form factor. So, below we compute this form factor using the 
standard approach of perturbative QCD assuming the scalar meson is a two quark state. 

In order to fix the convention, we recall the form factor definitions : 
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TABLE III: Branching ratios of measured PS channel decays of B mesons 



^ we used Br{f^ — > 2tt) = 0.68 ,3i in order to get the Br(B '^^ K •''/'^) from pubhshed results. 
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(M2(p2)|i"|Mi(pi)) = ipi+P2 



pAhM 



with q = Pi — P2- Projecting the amphtude on q one obtains: 



quiMMlL^lMM) = iml-ml)F,''^^''\q^) 



qu{MMMM\L-\Q) - {mi-ml)F^"^^'''^{q^) 



(11) 
(12) 



(13) 
(14) 



PQCD contributions to both amphtudes have exactly the same structure, so we compute qu{M2{p2)\L^\Mi(pi)) 
following ref. 26] then qu{^l2{p2)\L'^\Mi{pi)) is obtained just changing the sign oip2- The form factors are expressed 
in terms of the distribution amplitudes: 



^^{x,p) = 

'^ao{x,p) = 

where Iq is the identity in color space, p — ^^p^ and 



ilr. 



Ic 



4i^{x){p + m^)75 

{x){p + TOao)75 



(15) 
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4>Tr{x)dx 

(l)ao{x)dx 
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2^/27v: 



(17) 
(18) 



The wave functions 4'-K,ao {x) are given by |2 



4>T,{x) 



2N, 



:Ux{l - H 

/aox(l-x) (l + BiCl^\2x-l 



(19) 
(20) 



where = ISOMeV, fao — ^MeV,\Bifao\ — TbMeV, and C'l^^(2x— 1) is the Gegenbauer polynomial. Thus, the 
matrix element is expressed as: 

q^{nip2)\L^\aoipi)) = -C{R)^^gl / dxdy4>aMMy) 



Tr 



75 (P2 + "i^)7^Pi/g^i''(pi + maoh" 



Tr 



75 (P2 + ■rn^)qtiL'^' P21Y [pi + 'mao)^u 



(21) 



where C{R) = 4/3, [pi - P2Y = q^ ^ rn%, k = -xpi + (1 - y)p2, Pu = k + yp2, P21 = -fc + (1 - x)pi and 



752 22i2i/'2 2 2 \ 

P21 = i'^-y)m% + yml^-mly{l-y). 



(22) 
(23) 
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FIG. 1: Allowed region for F(f " and F^"" at one cr using experimental data on Br{B'^ n a^) and Br{B 



7T-a°) 



Integrating numerically, one gets 



^0 yj^B) 



0.004 



(24) 



It is important to notice that the CP-conserving phase of the annihilation contributions is not fixed since we only 
known the absolute value of Bi . 



GENERAL FRAMEWORK TO PREDICT B ^ PS 

One can proceed along similar lines for processes involving oq or /q scalar mesons. Instead, we use isospin, SU (2) 
quark symmetry and the quark contents of the scalar mesons to obtain relation between the form factors. We also 
used available experimental data to obtain constraints on the form factor values. It turns out that the consistency of 
the two sets of values so obtained provide further confidence on the approach. 

We assume the conventional quark content of the pseudo scalar mesons |^ and parameterize the mixing in the scalar 
sector in the strange-nonstrange basis as: 



cr = cos (psnn — sin (j)sss, 
fo = siiKpsnn + cos (j>sss, 



(25) 
(26) 



where nn = {uu + dd)/^/2, and the singlet-octet mixing angle 6s is related to (t>s by (t>s — Os — cos^^[l/\/3] ~ 55°. A 
diagrammatic analysis of the contributions of the form factors, which involve the quark composition and isospin and 
SU(2) symmetry between up and down quarks, lead the following relations : 
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(27) 



From these relations it follows that IF^^"] < IF^'^^I/V^. 



7 




It is possible to obtain similar relations between the annihilation form factors (Fq"'" and Fq ), however we will 
not use SU{3) symmetry since large deviations from the symmetry limit are expected. Using the experimental results 
given in Table Ijllip . the values for the scalar masses given in ref.j^l and the values given in table for the other 
form factors appearing in the amplitudes, it is possible to determine separately for each scalars uq and /o the allowed 
regions for the values of the form factors F^°^ {F^"^) and i^^f "° (-^cf ^) respectively. 

The results are summarized in figures Q and [3 Assuming that Perturbative QCD leads us the right order of 
magnitude for F^""^ , it follows that: 

0.14 < iFo^""! < 0.21 (28) 

We should note that this result is compatible (even if slightly smaller) with the predictions for F^ ° (0) — 0.55 ±0.22 
[iif obtained in a model-dependent way. Using Ea. H27|) . it follows that 

|-F(f°-^°| < 0.20 (29) 

Figure (|2Jl shows the values allowed by the experimental data when one standard deviation is considered. One observes 
that \Fq ■'"l < 0.20 requires a large contribution from |Fq |. In fact, the smallest value for |Fg |is around 0.05, 
which is more than one order of magnitude bigger than the PQCD evaluation of jF,^"'^!. It is interesting to note in 
this respect that: 



IF, 



foK, 



I _ '"-K 



\Fr\ ml 



12 (30) 



Scalar mesons as gqqq states. 



Several models where the scalars are four quark states jl3l l29l I3C 
at present time. We shall apply our method to one example, following 
scalars is given by 




have been published but no model is favored 
we assume that the quarks contents of the 
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FIG. 3: Branching ratio for Br{B^ ^^o^o )■ dot-dashed hne correspond to the four quark assignement for ao, and the 
space between both solid horizontal line is the value for Br{B'^' — > 7r"'"ag ) expected from 2 quark models for ao 



— uuds, Oq = dsus, aj] — —= (usus — dsds) , Kq — udds, Kq — udus, Kq = usud, Kq — dsud 

v2 



/o 



cos ( 



- -|- sdsd}j + sin udud, a = — 



sm ( 



(susM + sdsd) + cos udud 



(31) 



V2 ' ' " ' ^/2 

■where the mixing angle is obtained from the relation tan;/) = —0.19 (for 771^ = 0.47 GeV), so = —5.4° and 84.6°. 

It is "well-kno'wn that perturbative QCD predicts that the form factor will go like l/g^*^"^^) where n is the number 

of constituents of the hadron. If n = 4, (to^) is strongly suppressed and annihilation can be neglected. Varying 

the experimental results within one a and using the fact that in four quark models for scalars the annihilation does 
not contribute to the processes {F^"'^ = F^^° = 0), one concludes that: 



0.70 < F^^" < 0.75 



Proceeding in the same way with processes B a^Tr , o^tt^, one gets 



0.15 < IF, 



Baa I 



< 0.20 



(32) 



(33) 



which are closed to the values of l^jf""! obtained assuming the scalars are two quark states. 



Sub-dominant processes Br{B'^ t^^o.q ) and annihilation 



Once the allowed values for the form factors Fr 



Bfa 



and F^°'° have been constrained, we turn our attention to the 



subdominant processes Br{B^ tt^Gq) which is strongly suppressed by G parity and isospin. In Q the author 
concludes that a positive identification of this process is an evidence against the four-quark assignment of ao or else 
for breakdown of perturbative QCD. Using our estimates for annihilation contributions obtained using PQCD one 
can predict the values of Br{B^ Tr+aQ^). The results are presented in figure O In the four quark models for oq 
where annihilation is strongly suppressed, one gets for Br{B'^ tt+Oq ) « 10^^. In the two quark model for oq, the 
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main source of uncertainty is the phase of the annihilation contributions which cannot be fixed using PQCD. Varying 
between and tt the CP-conserving phase for the annihilation, it follows that: 

10"^ < Br{B° n+a^) < 4 x 10"^ (34) 

where the lower limit is obtained when annihilation interferes destructively and the upper limit when annihilation 
interferes constructively with the other contributions. Our conclusion is that this B decay cannot be used to distin- 
guish between two and four quark assignment of the ao, unless one can obtain an independent determination of the 
annihilation phase. 

Another channel suppressed by G parity is ^ Tr'^ao". It is interesting that this channel is better to distinguish 
between the 2 or 4 quark models for cq. Indeed, in 4 quark models for ao, with the value for iFff""] determined in 
previous sections, one gets 

2 X 10"^ < Br{B- 7r°ao ) < 10"^ (35) 

while in the 2 quark model: 

6.4 X 10"* < Br{B- 7r"ao ) < 2.4 x 10"'^ (36) 



CONCLUSIONS 



In this paper we consider processes for which the dominant contribution is suppressed. Using the factorization 
approximation and available experimental data, we estimate the effect of annihilation penguins contribution to the 
processes Br{B^ — n^aj) and Br{B^'~ K'^'~fo). We show that a consistent picture can be obtained when the 
scalars are described as two quark states, although one requires an important contributions from annihilation penguins 
to i?r(J3°'" K^'^ fo). Within our analysis the four quark models for /o cannot be excluded. 

Applying our estimates of the annihilation contributions to suppressed processes like Br{B^ '^'^o-a) ^"^^ Br{B^ — s- 
tt+Oq'), we conclude that the positive identification of B^ '^'^'^o cannot be taken as evidence for the four quark 
assignment of oq. This is in contrast with ref. where the annihilation contribution is not quantified. Relevant for 
this conclusion is the ambiguity in the CP-conserving phase of the annihilation penguins contributions. Our best 
candidate process to distinguish the nature of oq scalar is Br{B^ '^^o.q ) where the predictions for 4 quark models 
are typically one order of magnitude smaller than 2 quark models. 

Using the mesons quark content, SU{2) quark symmetry and isospin we derive relations between the form factors 
Fq^" to F^"" for different charge states. One can extend the analysis to SU (3) however one expects large deviations 
from the symmetry limit. This restricts the applicability of our approach to the four quark states since in that kind 
of models scalars necessarily involve strange quarks. 
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APPENDIX 



Below wc list the X^^ expressed in terms of the form factors. We quote only those needed to compute the branching 
ratios given in the paper: 

^3o„+ = < n-\{du)L\0 X 4\{ub)L\B° >= U{ml - rnDpf"^ {ml) 

^IV = < «o I Wi|0 X 77+1(^6)^15° >= -U{ml - rnDFi'^^ml) 

= < %\{du)L\0 X A{ub)L\B- >= fa{ml - ml)F^~^\ml) 

X<^_ = <n\uu)L\0><%\{db)L\B->=^{ml-ml)F,^-^-{ml) 

Xb'so = <7r-\{ud)L\0><S"\{ub)L\B- >=U{ml-mlo)F^~'^\ml) 

Xiol- - < 5°p-|(du)i|0 X 0\{ub)L\B- >= -/s(m|„ - ml)F;^"^'{ml) 

= <a,TT^\idu)L\0><0\iub)L\B- >=~fBiml-ml)F;;"''\m%) 

Xros, = < K"S°\{sd)L\0 X 0\{db)L\B° ~fB{mlo - ml)Ff''"{ml) 

Xs~K- = < S°K-\{su)l\0 X 0\{ub)L\B- >= -fsiml - mj,)Ff''' (ml) 

^{a-.+)u = < %^^\i^^)L\0 X 0\idb)L\B° >= -fBiml- - ml)F^°''\ml) 

= < a+n-\{uu)L\0 X 0\{db)L\B' >= -fsiml- - ml)F^^''~ {m%) 

X§°so = <K°\{sd)L\0><S\db)L\B''>=fjio{ml-mlo)Fi''^''{ml) 
X^Igo = < K-\{us)l\0 X S0|(«6)l|B- >= /^(ml - m|o)Fo^"^°(m2) 

= <a-|J^|0><;r+MB°>=mJ,-^f^Ff-"(m^) 

Till) '''^u 

= < S°\dd\Q X TT-\db\B- >= ^^^£^(m| - m^)Fo^"-"(m|j 

? 9 

„ — _ ^ 771 — Tin n 

X^%o = <a^\du\Q><T:%b\B- >=mj -^—^F,^-- {ml) 

nil) iild 

Xf_^_ = <50|ss|0><i<:-|s6|B- >=mso/sc^^^^^^Fo^"^"(m|o) 

■° " 777,5 — rri's 

Xfofio = < S\ss\0 X K'\sb\B' >= mso/so "'^^"'^> o^"^°(m|o) = (37) 
where Sq is a neutral scalar (og or /°). 



FORM FACTORS DEFINITIONS AND CONVENTIONS. 



In order to compute the amplitude using the factorization, we use the following parametrization of the form factors. 
The decay constants are defined as: 
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{0\A^\P{q)) = tfpq^. (38) 
(Ok-i7592|P(g)> ^ Zlka^L ^ f-pmp (39) 

mi + 7712 

(ao |J7^u|0) = /aoPp (40) 
(a(7|dM|0) = maofaO (41) 

Using the equations of motion (—j9''(9i 7^x7592) = ("t-i + m-2)Qi75'?2 and —id^{qi^^q2) = (mi — mT^f^i^jfT jlslEnj l on 
can show that fs = msfs/ (^^ii — ^2) and that /50 = for a neutral scalar. 
Form factors are defined as: 

< MM\L,\MM > = (pi +P2 - ^^^^-^q) ^^^^'^ + ^^^^^q,F,''^''%q') (42) 
<M2(p2)Mi(pi)|L^|0> = (p2-pi- 



9 ^ ; + ' 9 
9 / ^ 9 



/ 2 2 \ 22 

/ 1712 ~ nil \ EiA/aMi/ 2\ , m2 — mi j^iMaMw 2\ /^o^ 

(^P2-Pi ^2 lj ^+ ) + ^2 ^f^^o (9) (43) 

where L,, = 7^ ^-9^ = "j'^Pl- A factor of — i has to be added to the form factors in the case one of the mesons is 



— ' 2 

scalar. 



